Ferromagnetic instabilities in atomically-thin lithium and sodium wires 
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Using density functional theory the ground state structural, electronic, and magnetic properties of 
monatomic lithium and sodium chains with low average density are investigated. A metallic, zigzag 
ground state structure is predicted but, most interestingly, stable equilibria for chains under tension 
are predicted to be ferromagnetic, which can be traced to exchange effects arising from occupation 
of the second subband as a function of the interatomic distance. 
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Substantial variations in structural, electronic, and 
magnetic properties of solids have been predicted to oc- 
cur with decreasing dimensionality, and with the advent 
of new techniques in atomic manipulation and nanofab- 
rication [1-3], it is now possible to realize quasi-one 
dimensional systems experimentally. Using photoemis- 
sion spectroscopy, both Luttingcr-liquidlike properties [4] 
and Peierls instabilities [5] have already been observed 
in linear monatomic chains formed on insulating sub- 
strates. Despite confirmation of such unusual structural 
and electronic properties, the existence of magnetic or- 
der in quasi-one-dimensional systems remains somewhat 
controversial [6-8]. Generalizing an earlier prediction of 
Paul [9], the Lieb-Mattis theorem [10] formally prohibits 
ferromagnetic ground states for strictly one dimensional 
systems in the absence of spin- or momentum-dependent 
forces. Although strictly valid for single-band systems, 
this theorem is no longer applicable when two or more 
transverse subbands are occupied, and in fact we report 
here the results of first-principles calculations indicating 
a ferromagnetic ground state for lithium chains, as well 
as stable magnetic equilibria for quasi-one dimensional 
sodium chains under tension. The unexpected develop- 
ment of magnetic order in these classic nearly-free elec- 
tron systems can be traced to exchange effects arising 
from the changing occupation of a second subband as 
a function of the interatomic distance. Accordingly, at 
high [11,12] and now at low densities, lithium and sodium 
depart radically from the simple-metallic behavior that 
typifies their normal density electronic structure. 

To examine lithium and sodium chains over a wide 
range of interatomic distances, we use a plane wave im- 
plementation of density functional theory [13] within the 
local spin-density approximation [14], the Vienna ab ini- 
tio Simulation Package (VASP) [15] Our Li and Na ul- 
trasoft pseudopotentials [16] treat only the single outer s 
electron as valence; they also include non-linear core cor- 
rections [17], and both reproduce previously calculated 
structural and electronic properties of bulk Li and Na. 
We use a supercell containing four atoms aligned along 
only one direction (defining the chain) and include 20 
A of vacuum in the orthogonal directions to minimize 



interactions between neighboring chains. For several dif- 
ferent values of the chain-length per atom, x, we optimize 
both the interatomic distance, d, and the bond-angle, a; 
the resulting binding energies appear in Fig. 1. As can 
be seen, the energy has a pronounced minimum in both 
cases (around x = 1.7 A for sodium and x = 1.55 A 
for lithium), at which point d = 3.47 A and a = 61° 
for sodium and d = 2.85 A and a — 63° for lithium. 
Therefore we conclude that linear light- alkali chains are 
unstable to a lateral buckling into a planar zigzag or 
polymer-like configuration. Out-of-plane deformations 
do not result in lower energies. Previous experimental 
work by Whitman et al. [18] observed that monatomic 
cesium chains on GaAs and InAs substrates form a zig- 
zag structure, which then remained stable under tension. 
Recent first-principles calculations have also suggested a 
similar zigzag structure in monatomic gold chains [19]; 
a similar zigzag geometry was proposed [20] to explain 
expanded liquid Rb and Cs. 

Although both alkali chains satisfy the necessary con- 
dition for insulating behavior with two monovalent atoms 
per unit cell, band overlap nevertheless prevents the sys- 
tems from being insulating. When perfectly linear, the 
electronic structure of the each chain consists of a sin- 
gle half-filled s-like band, with two degenerate, unoccu- 
pied p-like transverse subbands lying above. The elec- 
tronic structure of the zigzag chain, however, exhibits 
two partially occupied subbands, one predominately s- 
like and one p-like. Formation of the zigzag structure 
lifts the two- fold degeneracy of the second subband (cor- 
responding to the two transverse p-like bands of the lin- 
ear chain), one of which then drops below the Fermi en- 
ergy, stabilizing the zigzag structure over the linear chain. 
This symmetry-breaking also increases the coordination 
of each atom from two to four nearest-neighbors, lowering 
Madelung-type terms and further stabilizing the zigzag 
arrangement. Additionally, we note that distances be- 
tween neighbors in both chains are longer than the bond 
length of the molecule (2.99 A for sodium and 2.71 A for 
lithium) but smaller than the shortest interatomic dis- 
tance in their equilibrium BCC solid phases (3.66 A and 
2.92 A respectively), reflecting an expected increase in 
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bond length with coordination number. 
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FIG. 1. Energy per atom (solid symbols) plotted as func- 
tion of the chain-length per atom, x, for (a) sodium and (b) 
lithium chains. The energy in each regime for both chains is 
well fit to i?(a;) = ao + a-ix -\- aiX'^ + a2,x^ + anx*' (solid, dashed 
and dotted lines). For each value of x, the ions are relaxed 
into their equilibrium positions until forces calculated from 
the Hellman-Feynman theorem are less than 1 me V/A. In all 
calculations we use a 30 x 3 x 3 k-mesh and plane wave cut-off 
of 20 Ry, together resulting in a total energies converged to 1 
meV/atom. When converging to the ground state, fractional 
occupancies are treated using a Gaussian smearing of cr = 0.1 
eV. The inset shows the three parameters characterizing the 
geometry of the chains: x (the chain- length per atom), d (the 
bond- length), and a (the bond-angle). 

Along with the energy minimum at a length of x = 1.7, 
the energy of the sodium chain has two noticeable kinks 
(at X = 1.9 and 2.2 A). Similar features are found in 
the lithium chain, though only one kink is observed (at 
X — 1.95 A) upon stretching (Fig. 1(b)). In both chains, 
we record complementary discontinuous changes in both 
d and a, each associated with major atomic rearrange- 
ments (as will be illustrated in Fig. 4 below). The addi- 
tional metastable equilibria observed in sodium exhibit 
surprising electronic and magnetic properties. As a; in- 
creases, the second subband becomes progressively less 
occupied, and this has significant consequences for the 
emerging magnetic properties, as we now show. (Al- 
though in the following our discussion will focus primar- 
ily on the sodium chains, our calculations indicate that 
equivalent conclusions are obtained for lithium.) The 
band structure corresponding to sodium chains for a; = 2 
A is shown in Fig. 2 and, remarkably, the bands for 
electrons with opposite spins are shifted with respect to 
one another, indicating that the sodium chains have a 
non-zero net spin polarization. 

Interestingly the lithium chain already has a net spin 
polarization at its equilibrium stabilizing length, thereby 
accounting for the absence of the first kink in Fig. 1(b). 
The spin-split p-like band, corresponding to antisymmet- 
ric orbitals oriented transverse to the direction of the 
chain, allows for electrons of a single spin character only 
and is thus solely responsible for the magnetic polariza- 
tion of the chain. The occupation number corresponding 
to this subband is so low (~ 0.15 e/atom) that the ex- 



change energy actually overcomes the excess of kinetic 
and correlation energies, resulting in the ferromagnetic 
state. We further note that the density of states (DOS) 
is observed to vary roughly as l/y^i? — (where is 
the energy of subband m with spin a at the F point), 
as is expected for free electrons in one dimension. This 
behavior arises from van Hove singularities at the bot- 
tom of each subband, which are especially prominent in 
one dimension. We note in passing that although our 
calculations carried out within the generalized gradient 
approximation (GGA) [21] result in a slightly larger in- 
teratomic distance d, they do not change the overall no- 
tably "polymeric" structure of the chains. The closing of 
the second subband and corresponding transition to the 
nonmagnetic state is found to occur at the same chain 
length, but the magnetic instability remains for a larger 
range of chain lengths within the GGA. 
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FIG. 2. (a) Band structure of the optimized sodium chain 
at a; = 2 A. Bands of opposite spin are denoted by the dot- 
ted lines. The Fermi energy (dashed horizontal line) is lo- 
cated at the bottom of the second spin-down subband. (b) 
Spin-resolved density of states (DOS). As expected for free 
electrons in one dimension, the contribution to the DOS of 
each subband follows E — eSi, being the energy of the 
bottom of subband m with spin o. 

The valence charge and spin polarization for sodium 
chain with x — 2.0 A appear together in Fig. 3. Notice 
that the charge is mainly confined to the region between 
the ions, and is almost uniform along the longitudinal 
direction of the chain. There does not appear to be any 
directional bonding between the ions: covalent bonding 
will therefore not account for the stability of the zigzag 
structure. The minima on either side of the charge den- 
sity maxima correspond to the location of the ions, where 
the valence charge density is small because of the pseu- 
dopotential approximation. Because of the p-symmetry 
of the second subband, the spin density also exhibits in- 
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teresting character: in Fig 3(b) two separate maxima can 
be seen on both sides with a minimum in the center. 
The spin polarization in this valley has opposite sign as 
compared to the net magnetic moment. The magnetic 
polarization of the chain is therefore highly anisotropic 
and, as previously indicated, the chain is ferromagnetic 
along X but has an antiferromagnetic character along the 
transverse direction of the chain. Similar charge and spin 
densities can be observed for lithium zigzag chains. 




FIG. 3. Valence charge density pup + Pdown (a) and spin 
density pup — Pdown (b) of the sodium zigzag chain at a; = 2 A 
as determined within the LSDA. The valence charge density 
maxima (red areas) extend along the chain between the ions, 
where the charge density is minimal (concentrated blue areas) 
owing to the pseudpotential approximation. In (b), oblate 
red regions of large spin density are located on both sides of 
the chain, away from the central charge density minima (in 
blue-green) . 

When the chain length is extended to x = 2.2 A , the 
second subband rises above the Fermi energy and the gap 
between spin-split subbands closes; the ferromagnetic or- 
der vanishes, resulting in the final kink observed in Fig. 
1(a) for the sodium chain. Wc may also note that for the 
paramagnetic state (having only one occupied subband) 
the energy is well fit to a polynomial function of x (as in 
Fig. 1), and the minimum energy obtained from this fit, 
corresponding to the cohesive energy of the chain with 
only the first subband occupied, is higher than the value 
calculated for the ground state of the chain. Thus filling 
the second band decreases the total energy and stabi- 
lizes the zigzag chain. The net magnetic moment per 
atom and the optimal geometry of the zigzag chains as 
a function of x for both alkali chains is shown in Fig. 4. 
The maximum value of the net magnetic moment of the 
chain (per conduction electron) is 0.18/iB ('^ 0.16/iB 
for lithium), and it decreases with increasing x; and for 
X > 2.2 A, only one subband is occupied and the chain is 
once again paramagnetic. The discontinuities in d, corre- 
sponding to the kinks observed for the energy and major 
atomic rearrangements, can be seen to reflect the onset 



of magnetic transitions. 

According to our results, the spontaneous magnetiza- 
tion in these zigzag chains originates from the exchange 
between electronic states in the second subband; these 
results are therefore not in violation of the Lieb-Mattis 
theorem, which holds strictly for systems in which only a 
single subband is occupied. Recent calculations [22] have 
shown that sodium quantum wires, described within the 
jellium model, undergo spontaneous magnetization for 
certain radii. In particular, they have predicted that the 
ground state is fully spin polarized when only the lowest 
subband is occupied. The effects of band structure and 
associated orbital hybridization, absent from the jellium 
model, clearly have an important influence on the stabil- 
ity of the chain. When the electron-ion interaction and 
atomic structure are explicitly considered, we have found 
that the ground state of the chain is paramagnetic when 
only the second subband is occupied. At larger chain 
lengths (i.e., x > 2.2 A ) the first subband is completely 
filled and although the Fermi wave vector makes contact 
with the edge of the Brillouin zone, we do not observe a 
Peierls transition. Instead, we find an intermediate state 
where the chain remains metallic and paramagnetic until 
X = 3.5 A , above which dimerization drives the metal- 
insulator transition. 
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FIG. 4. The magnetic moment per atom as a function 
of the length per atom, x, for both lithium (solid) and 
sodium (dashed) chains. Lithium chains exhibit a stable fer- 
romagnetic ground state; ferromagnetism develops in sodium 
chains only under tension. The inset shows the optimized 
bond-angle, a (solid line), and bond-length, d (dashed line), 
plotted as a function of the chain-length per atom, x. 

In summary, we observe that one-dimensional lithium 
and sodium chains are unstable to zigzag polymer-like 
structures. The transition to the zigzag structures in- 
creases the coordination number and lowers the energy of 
a minimally-occupied p-like subband (possessing a large 
van Hove peak). The spin-splitting of a p-like subband 
moves the van Hove peak to lower energies in lithium, and 
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a ferromagnetic ground state is predicted; for sodium, 
this ferromagnetic ordering is stable only when the chains 
are under tension. Our calculations clearly establish an 
intimate connection between the cxistance of magnetic 
order and zigzag arrangements. But the fact that mag- 
netic order is found at all in these structures raises the 
important possibility of the occurancc of more general 
ground states, such as, for example, non-coUinear mag- 
netic order, or even static spin waves as originally pro- 
posed for the alkalis by Overhauser [23]. Accordingly, 
future extension of this work along these lines would be 
of considerable interest. The results presented here for 
lithium and sodium may not necessarily be extended to 
chains composed of heavier elements. For the sake of 
comparison, monoatomic gold chains also show the same 
zigzag structures [19] and the energy shows two minima 
with the intermediate maximum associated to the clos- 
ing of a conduction band. However, we have found a very 
small magnetic instability in gold chains, /imax ^ 0.02/iB, 
and we believe that is related to the presence of d bands 
near the Fermi energy. The possibility of magnetic chain- 
or stripe-like structure is also predicted to arise in an im- 
portant system which might be seen to underlie the alkali 
chains considered here, namely the two-dimensional in- 
teracting electron gas at low densities [24] . Future experi- 
ments on inert substrates (e.g., cleaved rare-gas crystals), 
chosen so that they not interfere with the p-like subbands, 
would clearly be of very significant interest. 
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